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Abstract. The rheological behavior of soft glassy materials basically results 
from the interplay between shearing forces and an intrinsic slow dynamics. This 
competition can be described by a microscopic theory, which can be viewed as 
a nonequilibrium schematic mode-coupling theory. This statistical mechanics 
approach to rhcology results in a series of detailed theoretical predictions, some 
of which still awaiting for their experimental verification. We present new, 
preliminary, results about the description of yield stress, flow heterogeneities and 
activated processes within this theoretical framework. 



PACS numbers: 05.70.Ln, 64.70.Pf, 83.50.Gd 
1. Introduction 

The discussion of 'glassy' materials in textbooks is usually restricted to simple 
molecular glasses, such as silica or polymeric glasses, and hence covers the standard 
field of the 'glass transition'. An impressive number of experiments performed in 
the last decade shows that typical glassy effects are not restricted to structural 
glasses, but are observed in a much wider variety of experimental systems. As a 
result, 'glassy dynamics' is being actively studied in systems as different as dirty type 
II superconductors, disordered ferromagnets or ferroelectrics, disordered electronic 
systems, soft glassy materials, granular matter. These observations not only allow to 
draw interesting analogies between different systems, but also to use similar theoretical 
paths to describe various fields. This paper is concerned with the application of a 
theory initially developed in the context of the statistical mechanics of glasses to 
describe the rheology of soft glassy materials. 

The term 'soft glassy materials' was proposed in Rcf. [0 as a generic name 
for a large family of complex fluids which share a similar phenomenology: colloids, 
emulsions, pastes, clays... It was then assumed that their physical behavior basically 
results from the competition between an intrinsic glassiness, in the sense of large 
relaxation times, and shearing forces which 'accelerate' their dynamics. This 
competition was then theoretically described using the simple trap model of Ref. Q , 
phenomenologically extended to account for the effect of an external flow. The model 
was further studied in Refs. Since then, various phenomenological approaches have 
been proposed, which replace the concept of 'traps' by 'fluidity' Q, 'free volume' js], or 
'degree of jamming' but end up with extremely similar mathematical formulations. 
The interest of such models is that they are simple enough so that complicated flow 
and thermal histories can be easily implemented. The obvious drawback is that the 
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physics is put by hand from the beginning: shear is assumed to reduce the trap depth, 
or to increase the fluidity, the free volume, the degree of jamming. 

The approach we discuss here is less transparent but does not assume anything 
beyond the dynamical evolution of a Langevin type for a glassy system defined by 
its Hamiltonian 0|. Interestingly, it can alternatively be viewed as a nonequilibrium 
extension of schematic mode-coupling theories, as we describe in section ^. A very 
similar approach was recently followed in Ref. Q . Static and dynamic behaviors can be 
studied in such detail that a very complete physical description can be proposed. The 
theory is also sufficiently understood that its shortcomings, domain of validity, and, 
in principle, the path to possible improvements are also known. At present, only the 
case of steady shearing was investigated in great detail, as reviewed in section y. We 
hope to report on transient behaviors in a near future, the difficulties being mainly 
technical. The paper also contains two new steps. First, we discuss the issue of a 
yield stress, and present preliminary results for its behavior in section ^. Second, we 
show in section |^ that two dynamical solutions for a given shear stress are possible 
in a certain regime and discuss this feature in light of the recent observations of flow 
heterogeneities. 



2. A microscopic approach for nonlinear rheology 



What would be an ideal theory for glassy rheology? Ideally, one would like 
to start from a microscopic equation for the dynamic evolution of the system 
under study, say a supercooled liquid in a shear flow, and solve this dynamics 
exactly. This 'first-principles' approach is obviously an extraordinary challenge 
and approximations have to be made in order to get a system of closed dynamic 
equations. A well-known approximation in the field of the glass transition is the 
mode-coupling approximation Q which leads to the mode-coupling theory (MCT) of 
the glass transition [|lO[. That the resulting equations can be obtained in a standard 



perturbative way is discussed in Refs. ||9|, 11, |12|, although the original derivation 
makes this less transparent. 

Generally speaking, starting from an evolution equation for the density 
fluctuations dp{k,t), where fc is a wavector and t is time, and a given pair potential 
interaction V{k), the mode-coupling approximation amounts to a partial resummation 
of diagrams in the perturbative development of the dynamical action. One 
obtains dynamical equations (Dyson equations) which close on two-point correlation 
C{k,t,t') = {Sp{k,t)Sp{—k,t')) and conjugated response function At 
thermal equilibrium, the fluctuation-dissipation theorem Txik, t, t') = dC{k, t, t')/dt', 
implies coupled equations for the density correlators only: this is the MCT. The 
theory has then the pair potential (or alternatively the static structure factor 
S{k) ~ C{k,t,t)) as an input, and the dynamic behavior of the liquid, the 
correlators C{k, t, t'), as an output. 'Schematic' models focusing on a given 'important' 
wavevector, say ko, were formulated This amounts to write S{k) « l + A5{k — ko) 
and focus on C{t, t') = C(ko,t, t') as a single observable. These schematic formulations 
lose the 'ab-initio' character of the full MCT, but they do essentially capture the 
dynamical singularities arising when the wave- vector dependence is kept ]lO[ | . 

The MCT of the glass transition has been discussed at length in the literature |]l^, 
p^ . There are two points we would like to emphasize, though. 

(i) The 'perturbative' derivation described above where both correlation and 
response are kept, opens the door to the study of nonequilibrium behaviors. The 
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system can be out of equilibrium because the equilibration time scale is too large for 
the experimental window: one focuses then on aging behaviors . The system can 
also be driven out of equilibrium by some external force, for instance a shear force. 
This is the latter, rheological, situation we shall be interested in. 

(ii) Kirkpatrick and Thirumalai remarked that the schematic models can be 
exactly derived starting from some mean-field disordered Hamiltonian, like Potts and 
p-spin models p^ ]. The connection was studied further in an important series of 
papers 10 . Beyond the basic connection made between two fields (spin and structural 
glasses), these works more importantly complement MCT with all the knowledge one 
can get from the thermodynamic studies of, say, the p-spin model. It is worth recalling 
that the p-spin model exactly realizes most of the thermodynamic 'folklore' of the glass 
transition |p^ , such as such an entropy crisis, or the existence of many metastable 
states Both sides, MCT and mean- field disordered models, are now indissociable 
and form an ensemble that could generically be called the 'mean-field theory of the 
glass transition'. We shall see below how the knowledge of the free energy landscape 
(the 'spin glass part' of the theory) allows one to make qualitative predictions for the 
dynamics (the 'MCT part'), that appears to be verified when the actual calculations 
are made. 

3. Steady rheology: A brief review 

The microscopic approach to glassy rheology described in general terms in the previous 
section was quantitatively investigated in Ref. [Q. There, a schematic model (of the 
Fp family jl^) was extended to take into account the crucial ingredient that the 
dynamics is externally driven. Technically, this amounts to break detailed balance 
and leads to closed coupled equations for a correlator and its conjugated response 
function. As mentioned above, these equations can be alternatively derived from the 
driven Langevin dynamics of the p-spin model ]l8| . These dynamic equations were 
then solved in the plane (Driving force. Temperature) in the case of a constant driving 
force, and the results interpreted in the language of nonlinear rheology, the driving 
force being analogous to a shear stress a. Recall that in the absence of the shearing 
force, the model has a dynamic transition at a temperature Tc where the relaxation 
time diverges as a power law. Just above T^, one has the standard two-step decay of 
the correlation, characterized by functional forms described in detail in Ref. |l0|| . 

In the presence of the driving force, the whole physical behavior is encoded in the 
time decay of correlation and response functions. Their analysis leads to predictions 
at several levels [Q. At the macroscopic level first, one gets flow curves relating the 
viscosity r/ = J dtC{t) to the shear stress a, see Fig. |[ Beyond a linear regime 
restricted to the high-T, low-a region, the system is strongly shear-thinning. The 
power law ry cx 7"^/^ is obtained at T = Tc, where 7 = cr/77 is the shear rate. 
In the glassy phase, T < Tc, the model describes a power law fluid, 77 « 7^"'^-^^ 
with no 'dynamic' yield stress (a(T) < 1 lim^^o o'(7) = 0) Recall that the 

shear-thinning behavior is derived in the present framework without any assumption. 
The predicted macroscopic behavior compares very well with experiments and 
simulations |2^ . 

Second, in the spirit of MCT, predictions can be made concerning the functional 
form of the correlators in various time regimes. They are the so-called 'factorization 
property' at intermediate times and the 'superposition principle' at large times. Again, 
these are well verified in simulation works |23] . We are not aware of any experimental 




Figure 1. Flow curves ri = ri{a) for various temperatures above, at, and below 
Tc, in the model of Ref. M. Newtonian behavior is obtained at low a, high T. 
Power law shear-thinning is obtained everywhere else, with r/ oc ^"^/^ at T^. 



checks for the moment. 

Third, a modified fluctuation-dissipation relation between correlation and 
response is also derived. It reads 

where the effective temperature Tcs replaces the equilibrium temperature T p3| , |2^ . 
While the short processes show Toff = T, at long time one gets Tcs > T, showing that 
slow processes are 'quasi-equilibrated' at an effective temperature higher than the 
thermal bath temperature. Again, extensive numerical simulations have confirmed 
this prediction for various physical observables in a sheared fluid |2^, ^ |27[ |. 
Again, this important prediction has not been checked experimentally in this context, 
although specific protocols have been described in Refs. We refer to Ref. |^ 

for experimental studies of aging systems. 

4. Yield stress and activated processes 

Having access to a microscopic (spin glass) model behind the dynamical equations 
allows us to understand the evolution from the geometry of the corresponding 
phase space. At zero external drive, this connection has been studied in much 
detail 32 . Above the dynamical transition temperature Tc, the available 



phase space is dominated by one large basin in the free energy, corresponding to the 
'paramagnetic', or 'liquid', state. At Tc, a threshold level in free energy appears, below 
which the free energy surface is split into exponentially-many disconnected regions. 

The aging dynamics below Tc can be understood as a gradual descent to the 
threshold level, starting from high energy configurations |^^. The slowing down is 
the consequence of the decreasing connectivity of the visited landscape. See Ref. 
for recent similar statements at equilibrium. When the system is quenched from a 
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Figure 2. Existence of a static yielrj stress. The system is first prepared below 
the free energy threshold, as in Ref. |32j. The driving force is then applied. The 
two curves are both taken in the subsequent noncquilibrium steady state and show 
solid, (T < cry, or liquid, a > cry, behaviors; T = 0.1. 



high temperature, but at the same time driven by non-conservative forces, it remains 
similarly drifting above the free energy threshold, constantly receiving energy from 
the drive. 

On the other hand, if the system is prepared in one of the deep regions below the 
threshold, it remains trapped for all times From this picture, we expect in that 
case that a weak driving force will have essentially no effect beyond a trivial 'elastic' 
response of the system, as it is not strong enough to make the system overcome the 
barriers. If instead a strong drive is applied, the system should escape the low-lying 
valley and surfaces above the threshold, where the drive will suffice to keep it forever. 

We have investigated this situation using the method of Refs. |l8[ |34|. A typical 
result is shown in Fig. ||, which shows that the above expectations deduced from the 
topology of the free energy landscape are well verified when the calculations are done. 
This proves the existence of a static yield stress ayiT) in the model [l6[ |. Figure ^ 
is also clearly reminiscent of the thixotropic behaviors p5| ] commonly encountered in 
soft glassy materials. 

The difference between the free energy threshold and the equilibrium energy 
vanishes for T — > T~, and the very notion of a threshold disappears. One expects 
therefore, on physical grounds, that limj,^^,- aylT) ~ 0. This expectation is in 
qualitative disagreement with Ref. Q. We are currently studying the temperature 
dependence of aylT) in our model to conclude on this issue. 

We are now in a position to recall what is the main shortcoming of this mode- 
coupling type of approach. In any realistic system the structure of threshold and 
valleys may remain essentially the same, but now activated processes allow to jump 
barriers that are impenetrable at the perturbative level. In the spin language, barriers 
between states diverge at the thermodynamic limit, N ^ oo, where N is the total 
number of spins. 

Despite several attempts, the analytical inclusion of activation in this framework 
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Figure 3. Self-similar activated dynamics within the metastable states of the 
finite- Af version of the model driven by an external force. The first curve (top 
left) is taken in a single realization of the process. The 3 other curves are successive 
zooms of the first one, as can be seen from the time axis. A magnification factor 
of 6 is applied 3 times until elementary jumps are resolved in the last figure 
(bottom right). Parameters are TV = 20, T = 0.05. 



remains elusive. However, a qualitative understanding may be gained by studying the 
driven dynamics of the finite- TV version of the spin glass model corresponding to our 
theory. Keeping TV finite directly implies that barriers are finite, allowing thermal 
activation to play the role it cannot play in the thermodynamic limit |l^ . This type 
of study would not be possible within the MCT framework alone. We have thus 
investigated by means of Monte Carlo simulations the Ising p-spin model (for p = 3), 
defined by the Hamiltonian 

N 
i<j<k 

where Jijk is a random Gaussian variable of mean and variance VS/N. The driving 
force is implemented by the use of asymetric coupling constant, the amplitude of which 
being the driving force, which we call a, since it has the same role as the shear stress 
in Ref. Q; see also |2|]. We numerically find clear evidences for the existence 
of a critical driving force cry(T) below which the system is trapped ('sohd'), and 
above which it is not ('liquid'). The novelty is that activated processes now play an 
important role. This can be observed in Fig. || where the time dependence of the 
energy density is represented for a driving force amplitude just below the yield value. 
Strong fluctuations are observed in the energy density, and the system alternates via 
thermal activation between trapping periods and periods of freedom 18, p6[. 
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Figure 4. Distribution of trapping times ip(t) for various temperatures in the 
presence of activation. The different curves are vertically shifted for clarity. 
The full lines are the theoretical prediction, <fi{T) oc r~(-'^+^/^=' , with no fitting 
parameter. 



Interestingly, Fig. ^ also shows that this time evolution is self-similar, in the sense 
that zooming on a particular time window leads to a very similar picture. In Fig. ^ a 
zooming procedure by a factor of order 6 is repeated 3 times until unitary moves are 
resolved in the last figure. This strongly suggests that trapping times are power-law 
distributed. This is indeed what we find when the pictures of Fig. ^ are quantitatively 
analyzed to extract the distribution of trapping times (p{t). This is shown in Fig. ^, 
together with excellent comparison to the theoretical prediction (^(r) ~ r^Pi^) -^^rith 
an exponent P{T) = 1 -f T/Tc which is discussed below. Remark that the mean 
trapping time (t) — J dT'ip{T')T' is infinite for /3{T) < 2, which is equivalent to say 
that a < ay- Consistently with the previous discussion, we find that (jy{T) decreases 
when T is increased, and our preliminary results indicate that ayiT —>■ T~) — > 0. 
More quantitative studies are also in progress. 

This self-similar behavior is theoretically expected since a power law distribution 
can be obtained invoking activation dynamics in a landscape with exponentially 
distributed energy barriers, as is the case of the random energy model obtained 
in the p — > cxd limit. This view is obviously reminiscent of the dynamics of the 
trap model 0, the connection between the two approaches being mathematically 



described in Ref. |38 . These considerations in fact led us to predict the temperature 
dependence of the trapping times distribution in Fig. ^ where we took /3(T) = 
1 -I- T/Tc- Note also the similarity between Fig. ^ and recent numerical investigations 
of the 'potential energy landscape' at equilibrium around Tg, where results were 
quantitatively interpreted in the framework of the trap model |3S[] . 

5. Flow heterogeneities 

The results presented in the two previous sections have interesting implications that 
we now discuss in detail. We noted in section || that our model describes a power 
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Figure 5. A singular, or 'nonmonotonic', flow curve for T = 0.1 in the model 
of Ref. M. The horizontal dashed line is the value of the yield stress at this 
temperature. The region where spatial flow heterogeneities (shear-bands) may 
appear is delimited with the vertical line at 7 = 7c (see text). Vertical points 
symbolizes the fact that any shear stress a < cry is possible when 7 = (this 
value of 7 is of course not visible in a log-scale). 



law fluid at low temperatures, a ~ ^ the shear-thinning exponent satisfying 

I = a{Tc) < a{T) < 1. When the consideration on the yield stress of section |4| are 
included, the complete flow curve of the material admits a singular 7^0 limit, see 
Fig. ^. Hence, the flow curve becomes 'nonmonotonic'. This singularity is well-known 
at the experimental level Note that it physically results, in our theoretical 

description, from the existence of the threshold value in the free energy. Note also 
that a purely dynamical approach a la MCT misses this subtetly Q . 

A direct consequence is that there is a complete interval of shear stresses, 
cr G [0 : cry], for which two dynamical solutions are possible. This is examplified 
in Fig. 1^, where the correlation function is represented for the same value of the shear 
stress, but starting either from random initial conditions ('shear immediatly'), or from 
an equilibrated initial condition below the free energy threshold ('aging then shear'). 
Both curves are obtained in a driven steady state for the same value of the control 
parameters (T, a) , and are equally stable at the mean-field level. 

Looking again at Fig. ||, one sees that the flow curve defines a critical value 
7c of the shear rate, defined by (t(7c) = cry and represented by a vertical line in 
the figure. A very interesting question now is: what happens in an experiment if a 
global shear rate 7giobai < 7c is applied to the sample? An homogeous flow would 
indeed correspond to a stress cr < cry for which a second dynamical solution exists. 
The answer is known from experiments pO{ and simulations the system will 

spontaneously develop flow heterogeneities, where a flowing band (7iocai > 0) coexist 
with a non-flowing band (7iocai = 0), both supporting the same value of the shear 
stress. In that case, flowing regions will display a dynamic behavior described by the 
dashed line in Fig. ^ whereas the second, jammed, band will be described by the full 
line in Fig. o. This is also observed in a recent numerical simulation, see Fig. 2 of 
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Figure 6. Two possible dynamical solutions for T = 0.1 and cr = 0.1 in the model 
of Ref. [H. They coexist in a single sample when the system displays shear-bands. 
Compare with Fig. 2 of Ref. |4l[. 



Ref. |4l[] . This shear-banding phenomenon results then from the 'nonmonotonicity' 
of the flow curve in Fig. H, and can again ultimately be viewed, in our model, as an 
experimental consequence of the notion of threshold in the free energy landscape. 

An open problem, of course, is the dynamical selection of the bands in a real 
sheared material 



1 42 



43|. How does the system choose the relative size of the 
bands? The problem is difficult, since no thermodynamic argument (such as a Maxwell 
construction) can be applied to this nonequilibrium situation. The same question is 
presently much discussed in various complex fluids, like liquid crystals or wormlike 
micelles [Eal. 



6. Conclusion 

In these proceedings, preliminary results concerning the theoretical description of the 
yield stress, the flow heterogeneties and the role of activated processes in soft glassy 
materials via a nonequilibrium schematic mode-coupling theory were presented. An 
advantage of our approach is that no assumption are made, at odds with the more 
phenomenological models which are usually used in the field of rheology. This allows 
to make detailed microscopic predictions, much beyond the macroscopic rheological 
level where concurrent models are stuck. These predictions were briefly reviewed 
in section ^. We emphasize, as we did in Ref. |26 , that many of them are still 



experimentally unverified. In that sense, the situation is very similar to the mid-1980's, 
where schematic mode-coupling models were derived, but with little experimental 
confirmation of their main features. 

We mentioned at several places of this paper that our results were only 
preliminary, and we took advantage of this conference to present our 'work in progress'. 
Although more precise and complete studies are obviously necessary, we have shown 
that a qualitative understanding of the rheological behavior could be gained via this 
approach, beyond the steady rheological situation studied in Ref. @] . 
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A more quantitative approach would be to derive a 'full' nonequilibrium theory for 
sheared fluids. This program was recently undertaken by Fuchs and Gates ||]. Their 
derivation makes use of a projector formalism instead of the perturbative development 
described in section ^, and ends up with dynamic equations for correlators only. It is 
not clear at the moment if the qualitative discrepancies between the two approaches 
underlined in this paper are due to this alternative derivation where nonequilibrium 
effects (like the violation of the fluctuation-dissipation theorem) are not described. 
Comparisons between the two approaches would thus be very interesting. 

Last, we mention also that if flow heterogeneities are suggested in our model by 
the existence of two dynamical solutions for a given shear stress, this does not imply 
that shear-banding has to actually take place, nor does it allow to gain any insight 
into the problem of the selection between the two solutions. This certainly requires 
more ambitious approaches than ours. 
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